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In this work we investigate the X-ray edge singularity problem realized in noninteracting quantum 
dots. We analytically calculate the exponent of the singularity in the absorption spectrum near the 
threshold and extend known analytical results to the whole parameter regime. Additionally, we 
highlight the connections to work distributions and to the Loschmidt echo. 
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Introduction. In condensed matter theory the X-ray- 
edge singularity constitutes one of the most important 
paradigms appearing in a variety of different contextes. 
In the X-ray edge problem one probes the response of a 
fermionic system, interacting or noninteracting, subject 
to a sudden local perturbation. Its origin lies in the study 
of X-ray spectra of simple metals where it was shown that 
the absorption or emission of a photon corresponds to 
the sudden switch on or off of a local potential scatterer 
embedded in a noninteracting Fermi sea^ ^. Since then 
X-ray edge physics has been found in a variety of different 
systems such as Luttinger liquids with impurit}/^, Ander- 
son impurity and Kondo mo del j^^^, resonant tunneling 
current- volt age characteristics^^, gapped systems^ ^, de- 
coherence in two level systems^ or work distributions for 
weak local quenchea^. 

In this work we focus on the realization of the X-ray 
edge problem in noninteracting quantum dots allowing to 
study it in a controlled setup. The possibility to tune the 
system parameters in quantum dots enables to vary the 
relevant quantity in the X-ray edge problem, the phase 
shift of the conduction band electrons. We analytically 
calculate the absorption lineshape near the threshold of 
a suitably initialized quantum dot at zero temperature 
extending the known analytical results^ to the whole 
parameter regime. This is an important generalization of 
x-ray edge physics to an experimentally accessible setup 
and it constitutes one of the very few examples that al- 
low for exact solutions. We show that the absorption 
spectrum can be identified with a work distribution^^ for 
a local quench in a resonant level model. This identi- 
fication allows us to take a different point of view on 
the problem offering nice interpretations of the results. 
Moreover, we highlight the connection to the Loschmidt 
echo that can be related to the Fourier transform of the 
absorption spectrunP^EU. 

The paper is organized as follows. First, we introduce 
the experimental setup that allows to mimic the X-ray 
edge problem in quantum dots. Then we calculate the 
absorption spectrum near the threshold by an associated 
Riemann-Hilbert problem^^. In the end we show the re- 
sults and point out the relation to work distributions and 
the Loschmidt echo. 

Optical spectra of quantum dots. Below, we present a 
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FIG. 1: (color online) Schematic picture of a quantum dot 
coupled to a fermionic reservoir that allows to study x-ray 
edge physics in absorption spectra of quantum dots. The 
photon absorption of a suitably initialized quantum dot im- 
plements a sudden switch on of the tunnel coupling between 
dot and fermionic reservoir due to a capacitative coupling be- 
tween the excited electron and the residual hole. For details 
see text. 



possible experimental realization of X-ray edge physics 
in noninteracting quantum dots following the ideas of 
Helmes et al.^ and Tiireci et al.^. In Fig. [l]a schematic 
picture of the setup is shown. Consider a narrow quan- 
tum dot with a large splitting of the single-particle ener- 
gies coupled to an electronic reservoir. In order to decou- 
ple the two spin channels, a sufficiently large magnetic 
field h is applied to the dot, large enough to overcome 
the local Coulomb repulsion U between the spin-^ and 
spin-| electrons on the dot. Thus the problem separates 
into two independent electronic species and we consider 
without restriction a single channel of spinless fermions 
in the following. By varying the back gate voltage Vg 
the quantum dot can be tuned in such a way that the 
topmost occupied level lies far below the Fermi surface, 
{sf — £/i)/A ^ 1, and the lowermost unoccupied level 



far above it, {si — £f)/ ^ ^ 1, provided the level split- 
ting is large enough. Here, A = npoV'^ denotes the level 
broadening with po the density of states at the Fermi 
level and V the hopping amplitude of electrons between 
dot and reservoir. Thus, the lower level can be consid- 
ered as occupied and the upper one as unoccupied. If 
an incident laser beam with angular frequency uo excites 
the electron from the lower level into the upper one, a 
positively charged hole is left behind. We assume that 
the quantum dot is designed such that the local capaci- 
tative coupling [/eh between the excited electron and the 
hole is strong shifting the upper level £i to lower energies 
Sf to the vicinity of the Fermi level of the conduction 
band electrons, i.e., \£f — £f|/A < 1. The absorption of 
a photon creates a local exciton. The localized hole not 
only interacts with the dot electron, it also establishes 
a local potential for the conduction band electrons. Us- 
ing polarized laser light, it is possible to address just 
one of the two spin channels such that our description 
of a system of spinless fermions becomes directly rele- 
vant. Assuming that the hole is stable such that it can 
be considered as static, at least compared to the other 
time scales in the problem, we can model this system by 
the following initial (before absorption) and final (after 
absorption) Hamiltonians: 
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(1) 



The hole degree of freedom already has been integrated 
out and is contained in a constant energy shift /S.E of the 
final Hamiltonian. The operator c\ creates an electron 
with wave vector k in the reservoir. The colons : . . . : 
denote normal ordering with respect to the Fermi sea. 
For convenience, we measure the single particle energies 
relative to the Fermi level, i.e., Sp = 0. The operator c^ 
creates an electron on the upper level of the quantum dot 
whose energy differs depending on if a photon has been 
absorbed or not due to the local exciton energy. 

The Hamiltonian in Eq. (IT]) without the potential scat- 
tering term was introduced in the context of the X-ray 
edge problem by Kotani and Toyozawa^ ^^ to describe 
the X-ray spectra of metals with incomplete shells. They 
solved the problem analytically in the vicinity of the 
threshold for the case where the final local level lies far 
above or below the Fermi energy. A similar problem at fi- 
nite temperatures has been investigated in the context of 
decoherence in charge qubits^^. The combined influence 
of a local potential scatterer and a virtual bound state 
was first discussed by Kit a et al.^ who solved the problem 
analytically for the case where the final local level energy 
lies above the Fermi level, i.e., Sf > 0. 

The aim of this work is to extend the known ana- 
lytical zero temperature results to the whole parameter 



regime with a general framework that can also be useful 
in other contexts. This includes, for example, decoher- 
ence in charge qubits coupled to a defect leveP^'^ or 
absorption lineshapes of the Kondo exciton^^. 

Absorption spectrum. Assuming that the coupling be- 
tween the system and the light field is small, one ob- 
tains for the absorption spectrum A{uj), the rate at which 
photons are absorbed, in second order of the coupling 
(Fermi's golden rule) at zero temperature 



A{uj) = i^^Uen 
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Here, I'^o) denotes the ground state of the initial Hamilto- 
nian with energy Cgs and \en) is a complete orthonormal 
eigenbasis of the final Hamiltonian with corresponding 
energies e^. The constant prefactor K: contains the ex- 
perimental details such as the intensity of the incident 
laser beam and the system-light field coupling. Repre- 
senting the (5-function by an integral over phase factors 
one can relate A{uj) to a dynamical correlation function 
G{t) via Fourier transformation 



A{u) 



dt 

2^ 



J{uJ-£i)t 



G{t) 



(3) 



with 



G{t) = {0\e 



iHit^-iHft 



|0). 



(4) 

Here, |0) = c^lV^o) denotes a product state of the Fermi 
sea for the conduction band electrons with a filled lo- 
cal d orbital. In view of the X-ray edge problem, G{t) 
is the equivalent to the core-hole Green's function. The 
dynamical correlation function G{t) in Eq. Q is an im- 
portant quantity also in other physical contexts. The 
quantity C{t) = \G{t)\'^ is the Loschmidt ec ho tha t al- 
lows to quantify the irreversability of a systenP^EU, here 
Hi, under a perturbation, here Hf — Hi. Moreover, G{t) 
is the characteristic funtion of a work distribution P{uj) 
for a quench from Hi to Hf where P{oj) = hi~^A{uj) is 
the probability of having performed the work uj on the 
system under this protocoP^. The relation between ab- 
sorption spectrum and work distribution has been worked 
out recentl}/^^. This relation is evident from a physical 
point of view. The absorbed photon provides its energy 
u to the system. This is equivalent to having performed 
the work u. 

Analytic results for the dynamical correlation function 
G{t) in the asymptotic long-time limit t ^ oo have been 
obtained for the case where the final energy ej of the local 
d level lies above the Fermi level, i.e., Sf > G^. In the 
case without potential scatterer, Kotani and ToyozawsP 
calculated analytically the characteristic function G{t) in 
the limit where the final local energy level lies far above 
or below the Fermi level. In both systems, the long-time 
behavior of the dynamical correlation function G{t) is 
algebraic 



G{t) '^ {ir]ty 



(5) 



with an exponent 7 that only depends on the phase shift 
S of the conduction band electrons at the Fermi level in 
presence of the local perturbation. The prefactor rj of 
dimension energy is a high-energy scale of the order of 
the bandwidth. Due to the Friedel sum rule, d/n is the 
screening charge that determines the exponent according 
to the rule of Hopfield^^. 

In the following, we will extend the known results to 
the whole parameter regime including also the case where 
£f < 0. Although the problem is in principle quadratic, 
the mathematical difficulty stems from the fact that in 
contrast to the original X-ray edge problem an additional 
dynamical degree of freedom, the local d level, is switched 
on. Due to the quadratic nature of the problem, the final 
and initial Hamiltonians are both bilinear in fermionic 
operators, the characteristic function G{t)^ that is a ther- 
mal expectation value of exponentials in Hi and Hf, can 
be reduced to a single-particle problem. As has been 
shown recently, functions such as G{t) can be represented 
as determinant j^^ * ^^ -^ 



G{t) = detM, M = 1 - / + /^, 



(6) 



of matrices in the single-particle space. The matrix S 
with matrix elements 

Sw = {\ciScl\). S = e'''^'e-'''f\ 1,1' = k,d, (7) 

where |) is the true vacuum without any fermion, is es- 
sentially determined by the single-particle subspace of 
the S-matrix S. Their matrix elements reduce to the re- 
tarded Green's functions of the final Hamiltonian up to 
a phase. The initial state is encoded in the matrix /: 

fdd = 1, fdk = fkd = 0, fkk' = ^kk'0{ — k). (8) 

It will be convenient to separate the dynamics of the 
additional dynamical degree of freedom, the local d level, 
from the dynamics of the conduction band electrons. For 
that purpose, we write the matrix M in a block notation 



M 



A B 
C D 



(9) 



where 



A = Mm. Bu = Mdk. Ck = Mm. Duw = M^w (10) 

such that one obtains by use of an elementary property 
of the determinant: 



G(t) = detM = {A- BD-^C) deiD. 



(11) 



The matrix D now only includes reservoir states such 
that deiD can be calculated with techniques known from 
the original X-ray edge problem. But the separation 
of the reservoir and d level degrees of freedom comes 
at the cost of finding the inverse D~^ of an infintely 
large matrix. Using a technique by d'Ambrumenil and 
Muzykantskii^^, however, the evaluation of the determi- 
nant of D is equivalent to finding its inverse D~^ . 



Auxiliary Riemann-Hilbert problem. The determinant 
of D as well as its inverse D~^ can be calculated by solv- 
ing an auxiliary Riemann-Hilbert problenP^. For that 
purpose, the determinant detD is first written as: 

In detD = Tr [\n{D) - / In 5] + Tr/ In S (12) 

where a term x^ — Tr/ln^ has been added and sub- 
tracted. The function x^ — i^Et contributes a term 
linear in t with dE the ground state energy difference 
between initial and final Hamiltonian. In the following, 
we will concentrate on the nontrivial contribution 

X^=Tr[\n{l- f + fS)-f\nS]. (13) 



Introducing a new set of matrices 

s{X) = exp(Aln6') 



(14) 



the logarithm appearing in the expression of x^ can be 
represented in terms of an integral over the new variable 
A: 



X 



I dATr [(1 -/ + /.)- V-/^-'] 
Jo 



ds 
d\' 



(15) 



Thus, the problem of calculating the determinant has 
been transformed into the problem of inverting the ma- 
trix 1 — / + /'5(A) where for A = 1 we have D = 
1 — / + fs{l). This inversion of an infinitely large ma- 
trix can be performed by solving an auxiliary Riemann- 
Hilbert problem. For this mapping onto the Riemann- 
Hilbert problem it is convient to pass over to a time rep- 
resentation where the Fermi-Dirac distribution, for ex- 
ample, takes the following formi^ 



/(T,rO 



1 
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(16) 



It remains to determine the time representation of the 
matrix s. As shown in Ref.^^, the Wigner transform 
5'(r, E) of Skk' can be related to the instantaneous scat- 
tering matrix provided the scattering matrix or equiv- 
alently the product p{E)T{E) of the density of states 
p{E) and the T-matrix T{E) are only weakly depending 
on energy. Then we have that 

S{r,E) = Jdep{e)SE^s/2,E-s/2e-''^ 

-^ l-e{r)e{t-r)i27rp{E)T{E) (17) 

is either the identity matrix for r ^ [0, t] or the scattering 
matrix for the conduction band electrons. For the long- 
time dynamics t ^ oowe expect that only the low-energy 
excitations near the Fermi level are important such that 
we ignore the detailed dependence of the S-matrix on 
energy and perform the following approximation for r G 
[0,t]: 



S{t,E) -^ S{T,E = 0) = e 



.2iS 



(18) 



Here, 6 is the phase shift of the conduction band elec- 
trons at the Fermi level. Thus, the approximated s{X) is 
diagonal in its time-representation 

5(r, rO = (jS{r -r'), a = e^'^\ (19) 

The matrix 1 — f-\-fs can be inverted by simply requiring 

(l-/ + /5)-i(l-/ + /s) = l. (20) 

In the time representation this equation becomes a sin- 
gular integral equation with a Cauchy Kernel due to the 
singular nature of the Fermi-Dirac distribution. Such sin- 
gular integral equations exhibit an exact solution via an 
auxiliary Riemann-Hilbert problem^^. Let Y{z)^ 2; G C, 
be an analytic function in the complex plane except along 
the cut r G (0, t) where Y is supposed to fulfill the bound- 
ary condition 



r+(r) = ar_(r). 



(21) 



Here, r+(r) = Y{t + iO) and r_(r) = Y{t - iO). The 
solution to this problem is unique provided Y{z) -^ 1 for 

z ^ 00: 



Y{z) = exp 



1 

2^ 



'd»tW 



(22) 



The knowledge of the funtion Y then allows to construct 
the solution of inversion problenP^ 

(1-/ + /s)-1(t,t') = 

= Sir - r') + Y^\r)f{r, r') [Y.{t') - Y+{r')]{23) 

such that one obtains the following formula for x^- 

27r Jo Jo d^ 



X 



-Y^{r)a-' 



da 
dX' 



(24) 



For the calculation of x^, it is necessary to introduce a 
regularization scheme where the integral over r is cut off 
near the boundaries of integration according to ^ z^q"^ 
and t ^ t -\- i^Q^ where ^0 is an energy of the order of 
the bandwidth. Using this regularization one obtains for 
the asymptotic long-time behavior of detD: 



detD *^ e''^' 



(im 



-{S/nf 



(25) 



From Eq. (23), we can deduce the inverted matrix D~^ = 
(1 — / + fs{l))~^ and thus the long-time limit of the 
generating function G{t) according to Eq. (pTj) 



G{t) *^ (irit)-'', 7 = (1 - S/nf. 



(26) 



in agreement with the known result for the case ef > 
Sp in Eq. ([5| and consistent with the Hopfield rule of 
thumb^^. Thus, the known asymptotic behavior extends 
to the whole parameter regi me as a lready shown in nu- 
merous numerical calculations^^^IIII] ^j^^g result consti- 
tutes one of the rare cases where it is possible to obtain 
exact analytical solutions. 
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FIG. 2: Absorption spectrum A(ijo) as a function of the inci- 
dent light frequency uj near the threshold frequency cJth for 
different final energies Sf of the quantum dot level. For sim- 
plicity we restrict to the case ^ = without potential scat- 
terer. Here, A = np^V'^ denotes the half width of the hy- 
bridized level in the quantum dot with po the noninteracting 
density of states at the Fermi level. The thick lines have been 
obtained by NRG calculations^. The thin lines show the an- 
alytic results that fit perfectly the exact NRG data in the 
asymptotic low frequency regime. 



Loschmidt echo. As already mentioned before, the 
characteristic function G{t) is also related to the 
Loschmidt echcP^^ 



C{t) = \G{t)f = 1(0 \e'^''e-'^f'\ 0>f 



(27) 



The Loschmidt echo quantifies the stability of motion in 
time of a system, in this case the Hamiltonian Hi^ under 
a perturbation Hf — Hi. Thus, for long times t Eq. (26) 



states that, no matter how small the local perturbation 
is, the time evolution of the state |0) with the final Hamil- 
tonian drives the system into a subspace of the Hilbert 
space that is orthogonal to the initial state. From the 
Anderson orthogonality catastrophe^ it is known that 
the ground state of the final Hamiltonian is contained in 
this subspace. The system as a whole, however, does not 
evolve into the ground state of the final Hamiltonian as 
their overlap is constant in time. 



Absorption lineshape. From Eq. (26), one can deduce 



the behavior of the absorption lineshape near the thresh- 
old analytically 



Aico) 



u;^u;th 



0{uj - ujth) (^ - ^th)^ 



(28) 



that shows the typical power-law singularity. The sin- 
gularity is a consequence of the singular behavior of the 
initial Fermi-Dirac distribution of the conduction band 
electrons at zero temperature. At finite temperatures, 
the singularity is cut off. In Fig. [2| NRG data for the 
absorption spectrum is shown. For light frequencies uj 



in the vicinity of the threshold, the analytical power-law 
results included as thin solid lines fit perfectly to the ex- 
act NRG results. In view of the equivalence to a work 
distribution, the existence of the threshold in the absorp- 
tion spectrum is evident. In the beginning, the system 
is prepared in the ground state of the initial Hamilto- 
nian. The minimum energy, i.e., work, that has to be 
provided to the system by switching on the coupling to 
the resonant level is the ground state energy difference 
between initial and final Hamiltonian. Thus, it is impos- 
sible for a photon of energy less than the ground state 
energy difference to be absorbed. The singular behavior 
of the absorption spectrum shows that the dominant ex- 
citations that are created by the absorption process are 
low-energy excitations in the vicinity of the Fermi level. 
Conclusions. In this work we have presented a general 
framework that allows to determine analytically the sin- 



gular threshold behavior of absorption spectra in quan- 
tum dots at zero temperature. This establishes an impor- 
tant generalization of x-ray edge physics to experimen- 
tally accessible environments that can be used to observe 
x-ray edge physics in a controlled setup. Moreover, we 
highlighted the correspondence of the spectra to work 
distributions and to the Loschmidt echo. The presented 
framework might also be useful in other contexts such as 
decoherence in charge qubits or the Kondo exciton. 

Acknowledgements. We acknowledge valuable discus- 
sions with J. von Delft. We thank M. Hani and A. Weich- 
selbaum for providing us the NRG data. This work was 
supported by SFB TR12 of the Deutsche Forschungsge- 
meinschaft (DFG), the Center for Nanoscience (CeNS) 
Munich, and the German Excellence Initiative via the 
Nanosystems Initiative Munich (NIM). 



^ G. D. Mahan, Phys. Rev. 163, 612 (1967). 

2 R Nozieres and C. T. Be Dominicis, Phys. Rev. 178, 1097 

(1969). 
^ K. Ohtaka and Y. Tanabe, Rev. Mod. Phys. 62, 929 (1990). 
^ A. O. Gogolin, Phys. Rev. Lett. 71, 2995 (1993). 
^ A. Kotani and Y. Toyozawa, J. Phys. Soc. Jpn. 35, 1073 

(1973). 
^ T. Kita, K. Ohtaka and Y. Tanabe, J. Phys. Soc. Jpn. 56, 

4609 (1987). 
^ I. Affleck and A. W. W. Ludwig, J. Phys. A: Math. Gen. 

27, 5375 (1994). 
^ R. W. Helmes, M. Sindel, L. Borda and J. von Belft, Phys. 

Rev. B 72, 125301 (2005). 
^ H. E. Tiireci et a/., arXiv:0907.3854 (2009). 
^° K. A. Matveev and A. I. Larkin, Phys. Rev. B 46, 15337 

(1992). 
^^ V. V. Mkhitaryan, E. G. Mishchenko, M. E. Raikh and 

L. I. Glazman, Phys. Rev. B 80, 205416 (2009). 
^^ D. Segal, D. R. Reichmann and A. J. Millis, Phys. Rev. B 

76, 195316 (2007). 
^^ A. Silva, Phys. Rev. Lett. 101, 120603 (2008). 
^^ P. Talkner, E. Lutz and P. Hanggi, Phys. Rev. E 75, 

050102(R) (2007). 
^^ A. Peres, Phys. Rev. A 30, 1610 (1984). 



^^ N. d'Ambrumenil and B. A. Muzykantskii, Phys. Rev. B 

71, 045326 (2005). 

^^ A. Kotani and Y. Toyozawa, J. Phys. Soc. Jpn. 35, 1082 

(1973). 
^^ A. Grishin, I. V. Yurkevich and I. V. Lerner, Phys. Rev. B 

72, 060509(R) (2005). 

^^ B. Abel and F. Marquardt, Phys. Rev. B 78, 201302(R) 

(2009). 
^° M. Heyl et a/., to be published. 
^^ M. Heyl and S. Kehrein, arXiv:1006.3522 (2010). 
^^ J. J. Hopfield, Comments Solid State Phys. 2, 40 (1969). 
^^ D. A. Abanin and L. S. Levitov, Phys. Rev. Lett. 93, 

126802 (2004). 
^^ I. Klich, in Quantum Noise in Mesoscopic Systems, ed. by 

Yu. V. Nazarov (Kluwer, Dordrecht, 2003); 
^^ L. N. Oliveira and J. W. Wilkins, Phys. Rev. B 32, 696 

(1985). 
^^ N. I. Muskhelishvili, Singular Integral Equations (P. No- 

ordhoff Ltd., Groningen, 1953). 
^'^ A. Weichselbaum and J. von Delft, Phys. Rev. Lett. 99, 

076402 (2007). 
^^ P. W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 



